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1 Introduction 



One of the challenges facing the Standard Model is an understanding of the strong CP 
problem [1]. The Lagrangian of quantum chromo-dynamics (QCD) contains a term 

Lqcd ^ Gfj.^G'^'", (1) 

where G'^'^ is the gluon field strength tensor, Qs the QCD gauge coupling, and 6 a dimension- 
less parameter. This term violates both P and CP invariance. In the full QCD Lagrangian 
involving massive quarks, 9 is not a physical observable, but 9 defined as 

9^9 + arg(detMu detMo), (2) 

is. Here Mu and M u arc the up-quark and the down-quark mass matrices. The observed CP 
violation in weak interactions indicates that the complex phase in the Cabibbo-Kobayashi- 
Maskawa matrix is of order unity. Without additional symmetries one would expect that 
the quark mass matrices, and hence their determinants, also contain large complex phases. 
Therefore one would naturally expect from Eq. (2) that 9 is order one. However, this is in 
gross violation of the experimental limits arising from the neutron electric dipole moment 
(EDM). Eq. (2) leads to a neutron EDM of order rf„ ^ 5x lO'^^ 9 c cm. Current experimental 
limit on dn is dn < 10^^^ e cm [2], which puts a constraint 9 < 10^^'^. Why 9 is so small, 
when its natural values if of order one, is the strong CP problem. 

An elegant and popular solution to the strong CP problem is provided by the Peccei- 
Quinn mechanims [3]. It is based on on a global axial symmetry, U{1)pq, which is sponta- 
neously broken. In this mechanism 9 is promoted to a dynamical field, 9 = a/ fa, with an 
effective potential for this field induced by non-perturbative QCD effects. Here a stands for 
the pseudo-Goldstone mode of the spontaneously broken U{1) symmetry, the axion [4], and 
fa the axion decay constant (equal to the VEV of the scalar that breaks the PQ symme- 
try). Minimizing the axion potential would fix the vacuum expectation value (VEV) of a, 
or equivalently 9 to zero, thus providing a natural solution to the strong CP problem. 

The axion decay constant fa is model-dependent. In the original Weinberg-Wilczek 
axion [4], the PQ symmetry is broken spontaneously at the electroweak scale. This model 
contains two Higgs doublets $i and $2 with VEVs vi and V2. Without the QCD instanton 
effects, spontaneous breaking of the global U{1) symmetry would imply that the axion is 
massless. Including non-perturbative QCD effects which break the global U{1) explicitly 
the axion would acquire a mass given by [4] 

~ ~ 100 keV (3) 

V 

where m,r and /tt are the pion mass and decay constant, and v — \Jv1-\-V2 ~ 246 GeV. 
The axion decay constant is fa ^ v in this model. Laboratory data, especially from the 
non-observation of the rare decay K ^ t: -\- a, excludes this simple weak scale realization of 
the axion [1, 2]. 

The laboratory bound on the axion coupling can be evaded by introducing a singlet scalar 
S into the theory carrying a non-zero PQ charge. When this scalar acquires a VEV, the 
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global U{1) symmetry will break spontaneously. The axion decay constant in this case is 
fa ~ (S), which can be much larger than the weak scale. With v, the couplings of the 
axion with the Standard Model fields will be highly suppressed, by a factor v/ fa, and the 
laboratory limits are evaded. This class of models are the "invisible axion models". If the 
spectrum of the theory contains the scalars S and two doublets $i and $2, we have Dine- 
Fischler-Srednicki-Zhitnitskii (DFSZ) axion [5]. If the scalar spectrum contains a single 
Higgs doublet and a singlet S but the fermionic spectrum is enlarged to contain a pair of 
quark-ant iquark fields, we have the Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion [6]. 
In the DFSZ axion model the SM fermions carry PQ charges and the axion potential is 
induced by the QCD chiral anomaly associated with the SM fermions. In the KSVZ axion 
model only the exotic vector-like quarks carry PQ charges and their QCD anomaly induces 
the axion potential. 

The PQ symmetry breaking scale is constrained by a combination of laboratory, astro- 
physical, and cosmological limits to be in the range 10^'^ GeV < /a < 10^^ GeV. If fa is much 
above the weak scale, axions produced in the interior of stars can escape freely, draining the 
star of its energy too rapidly. The limit < lO-*^*^ GeV arises from the stabihty of stars and 
supernovae [1]. If /„ exceeds about 10^^ GeV, the energy density in the coherent oscillations 
of the axions will over-close the universe, contradicting cosmological limits. The mass of the 
axion in these invisible models is extremely small, ~ f-Kni-,^/ fA ~ 10~^ eV. Axion with a 
mass in this range can constitute a good cold dark matter candidate [1]. 

Quantum gravitational effects can potentially violate the global PQ symmetry which is 
needed for the axion solution. These effects, associated with black holes, worm holes, etc 
are believed to violate all global symmetries, while they respect all gauge symmetries [7]. 
When a gauge symmetry is spontaneously broken, often a discrete remnant survives [8]. Such 
discrete gauge symmetries are left intact by quantum gravity. Any such quantum gravity 
violations of the PQ symmetry should be extremely small for axion to be a viable solution 
to the strong GP problem. For example, in the invisible axion models, the gauge symmetries 
would allow for a term in the scalar potential of the form 

en 

where S is the singlet field carrying PQ charge, Mpi is the Planck scale, and n is a positive 
integer. This term would induce a non-zero 6 given by 

^ ~ M A 4 

iWpi I^QCD 

where Kqcd — 0.2 GeV is the QGD scale. Since the neutron EDM sets the hmit 9 < 10~^°, 
using fa ~ 10^^ GeV one finds that n > 10 is necessary in Eq. (4) [9]. This is indeed a 
severe constraint on axion models. 

A possible way to avoid the problems associated with quantum gravity is to identify 
Peccei-Quinn symmetry as an approximate global symmetry associated with a gauge sym- 
metry. Attempts have been made in the past with some success to achieve this by extending 
the low energy (below the Planck scale) particle content of the invisible axion models [9]. 
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In this paper we show how it is possible to use discrete gauge symmetries to stabihze the 
axion solution without enlarging the low energy particle content. This is made possible by 
a discrete version of the Grcen-Schwarz anomaly cancellation mechanism [10]. Superstring 
theory when compactificd to four dimensions generically contains an anomalous U{1) sym- 
metry with its anomalies cancelled by a shift in the pseudoscalar partner of the dilaton field 
(the string axion field). This anomalous U{1) is broken near the string scale when some 
scalar fileds carrying U{1)a charges are shifted to cancel the Fayet-Iliopoulos term associ- 
ated with the U{1)a (so that supersymmetry is left unbroken near the string scale) [11]. 
Often a discrete Z^v subgroup of the U{1)a symmetry is left unbroken to low energies. This 
Zjv symmetry, being of a gauge origin, is protected against quantum gravitational violations. 

Prom an effective theory point of view the conditions for Green-Schwarz anomaly can- 
cellations are 

A A 

-^i ■'^gravity r 
h ^ ^ 

where Ai are the anomaly coefficients associated with G'j x U{1)a with Gi being one of the 
factors of the surviving gauge symmetry. 5gs is a constant that is not specified by the low 
energy theory alone, ki are the levels of the Kac-Moody algebra. For non-Abelian groups ki 
are (positive) integers, although for U{1) factors ki need not be integers. All other anomaly 
coefficients such as GiGjGk and [[/(l)yi]^xG'j should vanish. In the case of SU {3)cxSU {2)lX 
U{1)y X U{1)a being the surviving symmetry, the GS anomaly cancellation conditions are 

where A3 and A2 stand for the [SU{S)c]'^ x U{1)a and [SU{2)l]'^ x U{1)a anomaly coeffi- 
cients. There is an analogous expression involving the hypercharge Y, however since Y is 
not quantized, the associated level ki is not an integer in general and this condition is not 
very useful from an effective low energy theory point of view [12]. The [?7(1)^]^ anomaly 
can be cancelled by the GS mechanism, but this condition also has an arbitrariness from the 
normalization of C/(1)a- The anomaly in [[/(l)^]^ x U{1)y also does not give any useful low 
energy constraint. 

When the U{1)a symmetry breaks down to a Zn discrete symmetry, the effective low 
energy theory would satisfy the following discrete Green-Schwarz anomaly cancellation con- 
dition [12, 13, 14, 15]: 

^3 + mN/2 ^ A2 + m'N/2 

ks ~ k2 ^ ^ 

where m and m' are integers. Heavy particles which are chiral under U{1)a can acquire 
masses of order the U{1)a breaking scale through ZAr-invariant Yukawa couplings, the con- 
tributions proportional to the integers m and m' account for these heavy fermions. 

In what follows, we shall impose the conditions of Eq. (8), which should be sufficient for 
the surviving symmetry to be protected from quantum gravitational effects. 
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2 Stabilizing the DFSZ Axion 



The non-supersymmetric DFSZ axion model [5] introduces two Higgs doublets Hu and Hd 
and a Standard Model singlet scalar S. The Lagrangian of the model relevant for the 
discussion of axion physics is 

C = Qu'^Hu + Q<fHd + Le'^Hd + Lv'^H^ + z/V^ - \{H^HdS^ + h.c). (9) 

Here we have used a standard notation that easily generalizes to our supersymmetric ex- 
tension as well. Q, L are the left-handed quark and lepton doublets, u'^, d"^, e'^, u'^ are the 
quark and lepton singlets. The seesaw mechanism [16] for small neutrino masses has been 
incorporated in Eq. (9). 

Eq. (9) has three U{1) symmetries, as can be inferred by solving the six conditions 
imposed by Eq. (9) on nine parameters. These three U{1) symmetries can be identified as 
the SM hypercharge U{1)y, baryon number U{1)b and a Peccei-Quinn symmetry U{1)pq. 
If we denote the charges of {Q, u'^, d^) as (g, m, d), the symmetries can be realized as 
B — q — u — d, PQ — —d, Y/2 — q/6 — 2u/3 + d/3. The C/(l) charges of the various particles 
under these symmetries are listed in Table 1. 
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1/3 


-1/2 


1 





-1/2 


1/2 
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-1 


-1 




















PQ 








-1 





-1 





1 





-1/2 



Tabic 1: Y/2, B and PQ symmetries corresponding to hypercharge, baryon number and PQ charge 
respectively. The charges are assumed to be generation independent. 



After Hd, Hu and S fields develop VEVs, the global PQ symmetry is broken and the light 
spectrum contains a Goldstone boson, the axion. Non-perturbative QCD effects induce an 
axion mass [1] given by 

~ 0.6 X 10-^ eV 12!!^, (10) 

where fa ~ (5') is the axion decay constant. 

We now apply the Green-Schwarz mechanism for discrete anomaly cancellation to stabi- 
lize the axion from quantum gravity corrections. Even though the model under discussion is 
non-supersymmetric, the GS mechanism for anomaly cancellation should still be available, 
since SUSY breaking in supcrstring theory need not occur at the weak scale in principle. 
Since baryon number has no QCD anomaly, any of its subgroup will be insufficient to solve 
the strong CP problem via the PQ mechanism. On the other hand, the PQ symmetry does 
have a QCD anomaly, although with the charges listed in Table 1 it has no SU{2)l anomaly. 
Since hypercharge Y is anomaly free, we attempt to identify the anomalous U{1)a symmetry 
as a linear combination of PQ and B: 

U{1)a^PQ + iB. (11) 
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According to the Eq. (11) and the charge assignment presented in Table 1, we have for the 
anomaly coefficients for the U{1)a, 

As = [SU{3)]^xU{1)a^-1 

A, ^ [SU{2)rxU{lU = l7 ■ (12) 

If we identify 7 = — /c2/(3A;2), the anomalies in U{1)a will be cancelled by GS mechanism. 
Thus we have 

U{1)a^PQ-~B. (13) 

The simplest possibility is /c2 = /cs = 1, corresponding to the levels of Kac-Moody algebra 
being one. Normahzing the charge of the singlet field S to be an integer, Eq. (13) can be 
rewritten as 

U{1)a^6{PQ)-2{B). (14) 

The corresponding charge assignment is given in Table 2. As discussed earlier, since hyper- 
charge Y is anomaly free, one can add a constant multiple of Y/2 to the U{1)a charges, and 
still realize GS anomaly cancellation mechanism. The charges listed in Table 2 assumes the 
combination -|(6Pg - 2B + |y). As can be seen from Table 2, this choice of charges is 
compatible with SU{5) grand unification. 

Suppose that the U{1)a symmetry is broken near the string scale by the VEV of a scalar 
field which has a U{1)a charge of in a normalization where all U{1)a charges have been 
made integers. A symmetry will then be left unbroken to low scales. Two examples of 
such Z]\f symmetries are displayed in Table 2 for A^ = 11, 12. Invariance under these Zjv 
symmetries will not be spoiled by quantum gravity, it is this property that we use to stabilize 
the axion. 

Potentially dangerous terms that violate the U{1)pq symmetry of Eq. (9) are S'"'/Mpf ^, 
HuHdS*"^ / Mpy~'^ etc, for positive integers n, m. For the induced 9 to be less than 10~^°, the 
integers n, m must obey n > 10, m > 5. The choice of N — 11, 12 satisfy these constraints. 
Note that a Ziq discrete symmetry would have allowed a term S'^, which would be inconsistent 
with the limit on 6. Z^ symmetries with A^ larger than 12 can also provide consistent 
solutions. Since by construction, the U{1)a symmetry in Table 1 is anomaly-free by GS 
mechanism, any of its Zfq subgroup is also anomaly-free by the discrete GS mechanism, as 
can be checked directly. In the Zn model, for example, we have A^ — A2 — 4. Consistent 
with the Zii invariance, terms that violate the U{1)pq symmetry and give rise to an axion 
mass are S^^/Mp^, HuHdS*^/Mpi etc, all of which are quite harmless. We conclude that 
the DFSZ axion can be stabilized against potentially dangerous non-renormalizable terms 
arising from quantum gravitational effects in a simple way. 



3 Stabilizing the KSVZ Axion 

The scalar sector of the non-supersymmetric KSVZ axion model [6] contains the SM doublet 
and a singlet field S. All the SM fermions are assumed to have zero PQ charge under the 
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Table 2: The anomalous ?7(1) charge assignment for the DFSZ axion model. Also shown are the charges 
under two discrete subgroups Zn and Z12 which can stabilize the axion. 



global U{l)pQ symmetry. The Yukawa sector involving the SM fermions is thus unchanged. 
An exotic quark-antiquark pair \[' + ^ is introduced, which transform vectorially under the 
SM (so the magnitude of its mass term can be much larger than the electroweak scale), but 
has chiral transformations under U{1)pq. The QCD anomaly needed for the axion potential 
arises from these exotic quarks. The Lagrangian involving the singlet field and these vector 
quarks is given by 

AL = 5*^ + /i.c. (15) 

When S field develops a VEV, the PQ symmetry is spontaneously broken leading to the 
axion in the light spectrum. 

The global PQ U{1) symmetry is susceptible to unknown quantum gravity corrections. 
We shall attempt to stabilize the KSVZ axion by making use of discrete gauge symmetries 
with anomaly cancellation by the Green-Schwarz mechanism. The most dangerous non- 
renormalizable term in the scalar potential that can destabilize the axion is S'"/Mpf^, as 
in the case of the DFSZ axion. We seek a discrete gauge symmetry that would forbid such 
terms. 

In order for the Green-Schwarz mechanism for anomaly cancellation to be viable, the 
anomaly coefficients A2 and corresponding to the [SU{2)i]'^ x U{1)a and [SU{3)c]'^ x 
U{1)a should equal each other at the U{1) level. This would imply that the ^ + ^ fields 
can not all be singlets of SU{2)l. The simplest example we have found is the addition of a 
5 + 5 of SU (5) to the SM spectrum. Such a modification is clearly compatible with grand 
unification. The 5 contains a (3,1) and a (1,2) under SU{3)c x SU{2)l. We allow the 
following Yukawa coupling involving these fields: 

£DA555 + /i.c. (16) 

If we denote the PQ charges of 5 and 5 as and 0, invariance of Eq. (16) under a surviving 
discrete Z^v symmetry would imply 

+ + s=piV (17) 

where p is an integer. In this simple model, all the SM particles are assumed to be trivial 
under the PQ symmetry. The discrete anomaly coefficients arc then A3 = = |(0 + 0) = 
\{pN — s). Since A2 = A^, the gauge anomalies are cancelled by the GS mechanism. As 
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long as > 10, all dangerous couplings that would destabilize the axion through non- 
renormahzable terms will be sufficiently small. We see that the KSVZ axion can be made 
consistent in a simple way. 

We have also examined the possibility of stabilizing the axion by introducing only a single 
pair of fermions under the SM gauge group, rather than under the grand unified group. Let 
us consider a class of models with a pair of fermions transforming under SU{3)c x SU{2) l x 
U{1)y X U{1)a as 

*(3,n,y,V') + *(3,n,-y,V;) , (18) 

along with a scalar field 5'(1, 1, 0, s). The Lagrangian of this model contains a term '^^S as 
in Eq. (15) and its invariance under an unbroken Zj^ symmetry imposes the constraint 

ij; + j/; + S ^ pN (19) 

where p and are integers. Since the SM particles all have zero anomalous U{1) charge, 
the anomlay coefficients arise solely from the (^ + ^) fields. They are 

1 — Tl 

^3 = -{nip + ni>) ^ -{pN - s) 

^ (» - DMn + 1) ^ ^ (n - l)„(n + D _ 

The Green-Schwarz discrete anomaly cancellation condition, Eq. (8), implies 

2(-m + bm') 

s = pN+ , r — ^— (21) 

where h = k^/k2. 

By choosing specific values of the Kac-Moody levels, one can solve for s, the singlet 
charge. For instance, in the simple case when = k2 b = 1, 

2{m' — m) 



— 3n 



-N. (22) 



We have normalized all U{1)a charges to be integers, including s, so the unbroken Zjv 
symmetry will be transparent. 

When n = 2, ^' and ^ are SU{2) doublets. One can calculate the charge of S from Eq. 
(22) and determine the allowed discrete symmetries. For 6 = 1, the solution is s = mod N. 
This solution would imply that 5"" terms in the potential are allowed for any n, in conflict 
with the axion solution. A similar conclusion can be arrived at for h = 1/2. For other values 
of 6, the Zjv symmetry typically turns out to be too small to solve the strong CP problem. 
For example, if 6 = (2,3, 1/3,3/2), the allowed discrete symmetries are (^4,^7,^3,^5). A 
special case occurs when 6 = 2/3, in which case s is undetermined, since A^/k^ — A2/k2- If 
one chooses s > 10, the KSVZ axion can be stabihzed in this case. 

If the quarks \& and ^ are triplets of SU{2)l, stability of the DSVZ axion solution can 
be guaranteed in a simple way. For b = k^/k2 = (1,2,3,1/2,1/3,2/3,3/2), which are the 
allowed possibilities if we confine to Kac-Moody levels less than 3, we have from Eq. (21), 
the unbroken discrete symmetries to be {Zg, Z21, Z^^, Zq, Z3, Zi^, Z30) respectively. For all 
Zn with A^ > 10, the axion solution will be stable against quantum gravitational corrections. 
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4 SUSY extensions of invisible axions and a natural 
link between Msusy, fa and the /j, term 

The models of the previous sections can be easily generalized to incorporate supersymmetry. 
The DFSZ axion is a natural extension of SUSY, since supersymmetry requires the existence 
of two Higgs doublets. For the axion to be weakly coupled (or invisible), we also need a 
pair of singlet scalars S and S. As we shall see, there is a natural link between the SUSY 
breaking scale Msusy, the PQ symmetry breaking scale /„, and the supersymmetric //-term. 



4.1 SUSY DSFZ axion 

The superpotential of the model is taken to be 

W = Qu^u + Qd'Ha + Le'^Hd + Lu'^H^ + Mru'^u' + Ai ^"J'"^' + ^2 ■ (23) 

As discussed earlier, the VEVs of the singlet fields S and set /„ ~ 10^^ GeV. Consequently 
the Ai term in Eq. (23) will induce a //-term for the MSSM Higgsinos given by [18] 

/( = Ai-^ «i 10^ GeV. (24) 
Mpi 

Furthermore, we have the following relation for the SUSY breaking mass terms in terms of 
the axion decay constant: 

Msusy ^ t^- (25) 

Thus, three a priori unrelated quantities get linked in these models. 

To see the link between Msusy and /„, we need to minimize the scalar potential"^ which 
contains both the soft SUSY breaking terms and the F-terms: 



Here C is the soft SUSY breaking term corresponding to the non-renormshzable couphng 
in the superpotential of Eq. (23).^ Let us assume for simplicity that ms — ttis- Relaxing 
this assumption does not have any major effects. 

By minimizing the potential of Eq. (26) [17] we obtain 



2 C ± VC^ - 12ms^ 



■^As far as the VEVs of the S and S are concerned, we can ignore the HuHdS'^/Mpi term, which can 
affect the VEVs by only 10"" GeV. 

^Notice that C can be either positive or negative. We shall follow a notation where C is positive. By 
setting S = |<S'|e'" and S = \S\e''^, C would appear in the potential as 2A2C ^^^^ cos2(q! + /?). Minimizing 
this potential we obtain cos 2 (a + /3) = —1. 
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where = {\S\)'^ = i^S^ . This is the desired relation between /„ andMsusY since Eq. 
(27) impUes 



fa ~ {\S\) ~ V^PiMsusY ~ 10'' GeV. (28) 

The scalar spectrum of this model contains two Higgs bosons from S and 5, besides those 
in the MSSM, with masses given by 

2 -2A2a;2CMpi + 24A2 V ^ Ox^Azx^CMpi + 4OA2 V 
"^'^ ^ M^- ' "^^^ = Mp? (2^) 

where x is given in Eq. (27). Note that both scalars have masses of order Msusy- Their 
mixings with the doublet Higgs is very small, being suppressed by a factor v/x ~ 10^®. 

Among the pseudoscalars in S, one combination remains massless and is identified as 
the axion. (QCD anomaly induces a tiny mass for the axion.) The orthogonal combination 
has a mass of order Msusy, given by 



2 4A2CX 

"^^3 - (30) 



The fermionic partner of the axion, the axino (a), receives a mass after SUSY breaking 
from two sources [19]. The superpotcntial terms induce an axino mass of order VuVd/Mpi ~ 
10~^ eV. A second source of its mass is through a Lagrangian term 

jCd J d^e{S^Sf/Mpi^ — > aaF* {S*) /Mpi^ ~ m^usv/Mpi ~ 10"=^ eV. (31) 

The fermionic component of the heavy pseudoscalar 5*3 obtains a mass from the superptential, 
of order Msusy ^ 

J dF9{SSf/Mp^ — > SpSp (SS) /Mpi ~ MsvsySfSf ■ (32) 

The axino is then the lightest SUSY particle in this model. The lightest neutralino (the 
LSP in conventional SUSY models) will decay eventually into an axino. To estimate this 
decay lifetime, we first note that the doublet Higgsino-axino mixing angle 9 is of order 
Msusy/Mpi)^/^ ~ 10~^. Suppose that the lightest neutralino is mostly a Bino and that 
its mass is larger than the lightest Higgs boson mass. The decay B a + can then 
proceed with a rate estimated to be F ~ (3q;i/20)(^^/ tan/3)^M^. For reasonable values of 
the parameters we see that the lifetime is of order 10^^ seconds, so the decay will likely occur 
outside the detector. For some other values of the model parameters, the decay length could 
be of order cm, which should be measurable experimentally. 

Turning now to the stability of the SUSY DFSZ axion, as in the non-SUSY case we 
make use of a discrete Zjy symmetry. We start with the same U{1)a charge assignment as 
in Table 2. One difference in the SUSY extension'' is that there is an extra S singlet field. 

^The Lagrangian actually lias four U{1) symmetries, hypercharge, Baryon number, Peccei-Quinn and R 
symmetries. However, the U{l)u is broken by the gaugino mass term to a discrete subgroup. Conventionally, 
the gauginos are assigned R charge of 1, after the i^-component of a spurion field Z acquires a VEV, to 
generate gaugino masses, the R symmetry is reduced to a Z2. This Z2 will be related to a subgroup oi B — L 
symmetry [15]. 
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The charge assignment under a subgroup is displayed in Table 3 ior N — 22. The Zu 
symmetry of Table 2 has been embedded into a bigger group Z22 in order to prevent a bare 
mass term SS in the supcrpotential, while allowing a term (SS)"^. This enables us to relate 
the fi term, Msusy and /„. This Z22 symmetry forbids to high orders all dangerous terms of 
the type 

cm Qn—m 

, . (33) 

that can potentially destabilize the axion. 

A second point that distinguishes the SUSY DFSZ axion model from non-SUSY case is 
that the anomaly coefficients now receive contributions from the Higgsino and the gaugino. 
(The absence of SS term in W , with {SSY allowed implies that the Grassman variable 9 
must transform non-trivially under the Zj^ symmetry.) Note that the anomalies of Z22 are 
cancelled by Green-Schwarz mechanism. 

In Table 3 we have also displayed a Z2 symmetry that turns out to be an unbroken 
subgroup of B — L, which is clearly anomaly free. This Z2 symmetry serves as it!-parity 
needed for the stability of the proton in MSSM [15]. 





Q u" d" L e" Hd S S a 


{A2, A,) 


Z22 

Z2 


3 19 1 11 15 11 22 18 13 20 11 
1111112 2222 


(19, 8) 

(0, 0) 



Table 3: The charges of particle corresponding SUSY DSFZ model with a Z22 discrete gauge symmetry. 



4.2 SUSY KSVZ axion model and its stability 

In the SUSY version of the KSVZ axion model we assume the superpotnetial to be 

(SS)'^ 

W = Qu^Hu + Qd^Hd + Le^Hd + Lv^H^ + Mrv^v^ + + ^— ^ . (34) 

Mpi 

Here \1/ + are the new vectorlike fields with QCD anomaly. Wc use the same mechanism 
as in the DFSZ model to relate fa with Msusy and ^ terms. It requires that the mass term 
SS be absent in W , but the coupling {SSY be present. 

The simplest possibility is to assume that the feilds ^ + ^ belong to 5 + 5 of SU{5). 
Gauge coupling unification that seems to occur within MSSM can be maintained with the 
introduction of complete multiplets of SU{5). The charge assignment under a Zu discrete 
symmetry are shown in Table 4. As before, the anomalies of Z12 vanish via the discrete GS 
mechanism. The lowest order term that breaks the PQ symmetry is 5*^^, which is harmless 
for the axion solution. 

In Table 5 we have displayed a variant KSVZ axion model where the Green-Schwarz 
anomaly cancellation occurs at a higher Kac-Moody level with ks — 3, and k2 — 2. We have 
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Q u" L -^^^ S S a 


(>l2, ^3) 


Z2 


0000000 11 156 
1111112 2 2 222 


(5, 11) 
(0, 0) 



Table 4: SUSY KSVZ axion model realized from Z\2 x symmetry. 



extended the solution obtained for non-SUSY case to the SUSY case here. There is a single 
+ field, belonging to a (3, 3) representation. Exact i?-parity is also reahzed as a Z2 
symmetry, as shown in Table 5. 





Q L E^ Ea ^ S S a 


(^2, A,) 


^30 
^2 


29 1 14 15 

11111122 2 222 


(24, 57/2) 
(0, 0) 



Table 5: SUSY KSVZ axion realized at higher Kac-Moody level from Z30 x Z2 symmetry. 



In conclusion, we have shown how the axion solution of the strong CP problem can be 

stabilized against potential quantum gravitational corrections via discrete gauge symmetries. 
Green-Schwarz anomaly cancellation mechanism plays a crucial role in our models. Both the 
DFSZ axion model and the KSVZ axion model are realized in rather simple ways without 
enlarging the low energy particle content. In the supersymmetric extension of these models, 
we have found an interesting link between the SUSY breaking scale, the axion decay constant, 
and the SUSY-preserving //-term. 
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